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^ ■ Abstract 

We examine the large branching ratio for the process B — > rj'K from the standpoint 
^ ' of R parity violating supersymmetry. We have given all possible $ p contributions to 

B — > amplitudes. We find that only two pairs of A'-type I/Lp couplings can solve 
this problem after satisfying all other experimental bounds. We also analyze those modes 
where these couplings can appear, e.g., B^ — > ir^K , B^'° — > K*^ ,0 rf'> , B^ — > (pK^ etc., 
and predict their branching ratios. Further, one of these two pairs of couplings is found 
to lower the branching ratio of B — > (pK^, thereby allowing larger £ = This allows 
us to fit B ± -> uK± and J B ± — > ujtt , which could not be done in the SM framework. 



*Electronic address: 
' Electronic address: 
* Electronic address: 



debchou@mri . ernet . in 
b-dutta@rainbow.physics.tamu.edu 
akundu@mri. ernet. in 



I. Introduction 

Recently, the CLEO collaboration has reported the branching ratios (BR) of a number of 
charmless nonleptonic B — ► PP and B — > PV two-body decay modes where P and V denote, 
respectively, a pseudoscalar and a vector meson. Some of these modes have been observed for 
the first time and the upper bounds on the others have been improved [|],[2[]. 

Among the B — > PP modes, the branching ratio for B ± — > r)'K ± is found to be larger than 
that expected within the Standard Model (SM). This result has initated lots of investigations in 
the last one year ||-§|. This kind of unexplained puzzle also exists in the B — > PV modes where 
it is found that the branching ratios of B^ — > (pK^, 5 ± — > ujtt ± and B ± — > ujK^ are hard to fit 
simultaneously |||9|. Present attempts to explain the large branching ratio BR^B^ — > rj'K^) 
involve large form factors and/or large charm content for rj', with contribution arising from 
b — > sec — > sr] f (r]), and low strange quark mass [f§-|(|. In an interesting paper [0], consequences 
of large B — > rj'K branching ratio from purely SU(3) viewpoint has been studied. 

In this paper we try to address the large BR problem from the standpoint of -R-parity (R p ) 
violating supersymmetry (SUSY) theories. Motivations for invoking SUSY and its lfi p version 
have been discussed in detail in the literature [[ l]]. Some of its effects on I?-decays have also 
been investigated [0 . Since the new interactions modify the SM Hamiltonian, it is natural to 
revisit these calculations and try to see whether the above mentioned puzzles can be solved. 
We calculate the QCD-improved short-distance part with the usual operator product expansion 
and Wilson coefficients (WC), while the long-distance parts are calculated by the factorization 
technique which is very successful in estimating B — > D decays. The requirement that any 
"new physics" solution of the perceived anomaly does not overly affect other observables that 
are in good agreement with the SM predictions restricts us to two particular sets of couplings 
within the ]fl p scenario. Interestingly enough, we find that one of these sets also leads to a 
better fit for the decays B ± — > (f)K , B ± — > ujti^ and B^ — > ljK ± . 

We organize the letter as follows. In section II, we give a very brief introduction to the 
SM and Jji v Hamiltonian, and list the possible ]jl p operators that can contribute to charmless 
decays. We discuss the B —* PP and B — > PV decay modes in section III. The new physics 
contributions to the decay modes B ± — > rj'K ± and B^ — > rjK^ are shown. In section IV, we 
discuss how ljl p can raise the branching ratio of B — > rj'K without jeopardizing other decay 
modes. We make predictions about the yet-to-be-observed channels which can be tested in the 
upcoming B-factories. We also discuss how to fit the new results in B — > PV modes in presence 
of the new couplings which are used to raise the BR of B — ► rj'K. We conclude in section V. 
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2. Effective Hamiltonian for charmless decays 
2.1 SM Hamiltonian 

The effective Hamiltonian for charmless nonleptonic B decays can be written as 



G p 
V2 



12 



H = =^= V ub V: q J2 dO t - V tb V* J2 c i°i + h - c - (!) 



i=l,2 i=3 



The Wilson coefficients (WC), take care of the short- distance QCD corrections. We find all 
our expressions in terms of the effective WCs and refer the reader to the papers [[7|[l3|-|T5f for 



a detailed discussion^. We use the effective WCs for the processes b — > sqq' and b — > dqq' from 
ref. @. The regularization scale is taken to be fi = nib. In our subsequent discussion, we will 
neglect small effects of the electromagnetic moment operator Ou, but will take into account 
effects from the four-fermion operators 0\ — 0\q as well as the chromomagnetic operator On. 

2.2 The R p - violating Hamiltonian 

The superpotential of the minimal supersymmetric standard model (MSSM) can contain 
terms, apart from those obtained by a straightforward supersymmetrization of the SM potential, 
of the form 

W $p = KiLiH 2 + XijkLiLjEZ + KjkLiQM + K^Dpi (2) 

where E t , Ui and Di are respectively the i-th type of lepton, up-quark and down-quark singlet 
superfields, Lj and Qi are the SU(2)^ doublet lepton and quark superfields, and H 2 is the Higgs 
doublet with the appropriate hypercharge. Symmetry properties dictate that = —Xjik 
and \'lj k = —X'ikj- Apparently, the bilinear term can be rotated away with a redefinition of 
lepton and Higgs superfields, but the effect reappears as As, A's and lepton-number violating 
soft terms ]TB| . The first three terms of eq.(0) violate lepton number whereas the fourth term 



violates baryon number. Thus, simultaneous presence of both sets would lead to catastrophic 
rates for proton decay, and hence it is tempting to invoke a discrete symmetry which forbids 
all such terms. One introduces the conserved quantum number 



R P = (-1 



,3B+L+2S 



which is +1 for the SM particles and —1 for their superpartners. However, to prevent proton 
decay, one needs to forbid only one set, and not necessarily both. This leaves us with the 
possibility of additional Yukawa interactions within the MSSM, many consequences of which 
have already been discussed extensively in the literature. 



1 Since the Ifip operators will be shown to be small, their mixing with the SM operators may safely be 
neglected at the current level of accuracy. 
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For our purpose, we will assume either A' or A"-type couplings to be present (A-type couplings 
do not lead to nonleptonic decays), but not both. Assuming all ljl v couplings to be real, the 
effective Hamiltonian for charmless nonleptonic 5-decay can be written as^ 



He/f{b —> djd k d n ) = df kn [d na ^djp dkpJnRb a ] + df kn [d na ^bp dfc/37/^j'a] 

He'ffib ~* UjU k d n ) = uf kn [u ko niUjl3 dnPlfiRba] , 

Heff(b ~> djd k d n ) = 2^jkn[dkalRdj/3 dnfilfiRba — dkcHlldja d n pJnR.bp\ > 

H e 'ff( b djdkdn) = u'j^UkaJ&Ujp d n/ 97/xfl&a ~ U ka ^U ja d n p^^ R bp] , 



(3a) 



with 



3 A' A' 3 A' A' 

,R \ -> A ijk A in3 iL \ A i3k A inj r • i _ -[ n\ 

d * kn ~ ^ 8m? ' d * kn ~ 8m? ' V'*'"- 1 ' 2 * 

1 = 1 ViL 1=1 ViL 

3 \i \i 

R \ ^ A ijn A ik3 

U jkn 



E^fV^> (;,* = !, n = 2) (3b) 



i=i 6m e lL 

3 \ll \/l 2 \ll \/l 
jll A ij3 A ikn /, v - ^ A ji3 A kin / ■ -\ o t, 1 o\ 

a jkn - 2-^~TZ5. ' U jkn - 2^~ZZ3. ' U - J-iA K - n - A )- 

i=l ^" L u iR i=l ^" L u iR 

where a and f3 are colour indices and 7^ L = 7^(1 ± 75). The parenthetical remarks on the 
subscripts concentrate on only the relevant couplings. 

As is obvious, data on low energy processes can be used to impose rather strict constraints on 



many of these couplings fl7H19|j • Most such bounds have been calculated under the assumption 
of there being only one non-zero ljl v coupling. From eq.fl3"a|), it is evident that such a supposition 
precludes any tree-level flavour-changing neutral currents, thus negating the very aim of this 
paper. However, there is no strong argument in support of only one fop coupling being nonzero. 



In fact, it might be argued [[Tj|] that a hierarchy of couplings may be naturally obtained on 
account of the mixings in either of the quark and squark sectors. In this paper we will take a 
more phenomenological approach and try to find out the values of such fo v couplings for which 
all available data are satisfied. An important role will be played by the type couplings, the 
constraints on which are relatively weak. 



3. B -> PP and PV modes 

We consider next the matrix elements of the various vector (V^) and axial vector (A^) quark 

currents between generic meson states. For the decay constants of a pseudoscalar (P) or a 

2 In this paper, we will not consider the CP-violating effects of these couplings, i.e., we will assume all of 
them to be real. However, the fact that they may not all be real leads to interesting consequences. 
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vector (V) meson defined through 



(0|A M |P(p)> 
(0K\V(e,p)) 



(4a) 



we use the following (all values in MeV) ||, 

f u = 195, f K . = 214, f p = 210, U = 134, f K = 158, = 1.10/,, f m = 1.34/.. (4b) 

The decay constants of the mass eigenstates rj and rj' are related to those for the weak eigenstates 
through the relations 

fs . n , fl n re n f% . n , /l 



f" = -i± s i n ^ + -i^cosfl 
" y/E y/3 

# = j=cos9-j=sin9, 



/- = _2^Lsinfl + ^cosfl 



/s /i 
2^= cos 9 = sin 9. 

y/E y/3 



The mixing angle can be inferred from the data on the 77 decay modes |20| to be 9 
Whereas the only nonzero B — > P matrix element can be parametrized as 



-22°. 



(P(p')\V,\B(p)} = [(p' + p),- 

the B — > V transition is given by 

(V(e,p')\(V,-A,)\B(p)} 
IV 



m 



B 



m, 



B^P 



+ 



m, 



ml 



B^P 



(5a) 



m B + niy 



+i 



[m B + m v )A 1 e* + e* • q { -A, 



m B + m v 



+ 2m v % (A - A 2 



(5b) 



^B^V 



with 2mvA 3 = (m^ + m v )Ai — (m B — m v )A 2 . All of the quantities Fq , y and Fq { 
have a formfactor behaviour in q 2 = (p — p') 2 - Note that F\ = F at q 2 = 0, and, to a very good 
approximation, we can set F(m 2 p 2 ) = F(0) for B decay formfactors since the q 2 dependence is 
dominated by meson poles at the scale m B . Flavour £77(3) then allows us to write 



^0,1 



B — >T]' 



sin 9 cos 9 



y/E \/3 



^0,1 


F 


P 0,l 


/ cosfl 


- F [VE - 



(5c) 



There seems to be considerable variation in the range of F estimated in the literature. Bauer 
et al estimate it to be 0.33 [21] while Deandrea et al get a value of 0.5 |[22|| . We find that while 
within the SM, the combination (77 = 0.36, [K^/Kftl — 0.07) yields a good fit to 77 — > tttt and 
B — > ttK data [[{J , introduction of 7# p interactions allows larger values of F. As for the 77 — > V 
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formfactors, it can easily be ascertained that, of the four, only A$ is relevant for the B — > PV 
decays that we are interested in. For the current W7 M (1 — 75)6, we have 

^0 = 7/^> ^0 = G, A = —j=, (5d) 

where we use G=0.28 ||. The only remaining parameters of interest is the mass of the strange 
quark for which we use m s (l GeV) = 165 MeV leading to m s (mb) = 118 MeV. 

3.1 B ± -> rt'irf)^ 

The effective SM Hamiltonian for this decay and its matrix elements are well-studied and 
can be found in Refs. 0,0. As for the ft p operators, it is easy to see that only six of them may 
contribute (with none from the u" set) and may be expressed in terms of 



Aau 



A 



(M 2 \Jg\B) (M 1 \J lll \0) 



where J and J' stand for quark currents and the subscripts b and / indicate whether the current 
involves a b quark or only the light quarks. Neglecting the annihilation diagrams*] we have, for 
the B — > 7]K matrix elements, 



M x ' 



(d 



R 
121 



u 112 



tA$ + (d 



L 

222 



d R 

u 222 



m . 
— (At 



AU 



+ (^121 ~~ ^m) ~r^v 



u 



R 
112 



2m 2 K A K 



(m s + m u )(m b - m u ) 



(6a) 



where m = mtjim^ — m s ), andfj 



M x " = d'Ul - £)A$. 



(6b) 



Analogous expressions hold for B — > rj K where we have to replace A™ by A™,, A^ by A s , 
and m v by m^. We note that \'[ 12 and A" 13 are bounded to be very small irrespective of the 
presence of other ft p operators, and hence may be neglected. For the numerical analysis, we 



take m 



rrirj and m Vl 



m v >. 



4. Analysis 

We are now ready to discuss our results. Our goal is to explain the branching ratio for 

the B^ — > rj'K^ decay while satisfying the experimental numbers (limits) for all other related 

decays (see Table [[]). To set the perspective, consider the solid curve in Fig. 0(a), wherein we 

3 Such processes cannot be treated under the factorization ansatz, but are expected to be negligibly small in 
any case. 

4 Note that (0\si l5 s\r)W) = _ % l0 )m* w /2m. |. 
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Mode 


5i? x 10 5 


SM theory xlO 5 


Mode 


5/? x 10 5 


SM theory xlO 5 


B+ _> V 'K + 
B + -> f]'K* + 
B+ _> ^+ 
5+ -> 77^*+ 


6.5l{j±0.9 

< 13 

< 1.4 

< 3.0 


0.8-4.3 
0.01 - 0.18 
0.06-0.14 
0.14-0.31 


B° -> r)'K° 

-w , Pi § -r . O 

B° -> r^ 
5° -> r?K *° 


4.7l|£±0.9 

< 3.9 

< 3.3 

< 3.0 


0.7-4.1 
0.03 - 0.18 
0.03-0.14 
0.1-0.5 




2 - 3 -i:o ± °- 4 

< 1.6 

< 2.0 


1.1 - 3.5 
1.0 - 1.4 
0.3 - 1.3 


-i— »n n T t- n 

B° -> 

nil 4- — 

5 U — > 7T + 7T 


< 4.1 

1.5lg;|±0.1 

< 1.5 


0.6 - 1.9 
1.1 -2.1 
0.8 - 1.5 


B+ -> 0/1 + 
5+ -> cuif + 


< 0.53 

i.5lg;I ±0.3 


0.07-5.0 
0.01-3.5 


5+ -> 


i.i!8:i ±0.2 


0.06-1.7 



Table 1: Branching ratios (or upper bounds) for various B -meson decays. Also shown are the 
theoretical predictions based on the SM only fftBjj . 

have plotted BR{B ± — > rfK ) as a function of £. It is quite apparent that only for very small £ 
could we hope to reconcile the SM predictions with the observations. One may argue, though, 
that such a conclusion is unwarranted in view of the uncertainty in other parameters such as 
F, the CKM elements V c b and V u b, the angle 7 of the unitarity triangle, and the strange quark 
mass 0. Consider instead the ratio BR(B^ — > rj'K )/BR(B — > ir + K°) which is independent 
of F and V c b- In Fig. |2], we plot this ratio as a function of 7 for £ = 0, so as to maximize it. 
Clearly, the SM prediction falls well below the experimental number (remember that 7 ~ is 
unable to account for the observed CP violation in If -system). 

Thus, if we demand that lji v solve the B —>■ r(K anomaly, the relevant operators need to add 
constructively to the SM amplitudes. We make a simplifying assumption here. Rather than 
consider the most general case, we restrict ourselves to exactly one non-zero product in eq. (paf) 
and discuss its consequences. This immediately restricts us to particular signs for each of the 
combinations. To wit, we need one of G?222> ^112 > ^ii2 ^^d d^ 2 to be positive. On the other hand, 
only negative values for the other four combinations d 222 , d± 21 , df 21 and d'{ l2 could explain the 
enhanced BR. We shall concentrate on only the first set. 

It is easy to see that itf 12 also enhances the BR(B — > n + K~). Since there exists a stringent 

experimental bound on this mode, the largest allowed value for uf 12 is too small to explain the 

BR(B — > rj'K). Similarly, the small enhancement due to d^ l2 , which occurs only for large value 

of £, is unable to explain the anomaly. Thus, we are left with only two terms, namely d 222 and 

5 The branching ratio of B — > r/'K increases slightly with the increase of the 77 — 17' mixing angle 9 0, but 
since the experimental constraint on this mixing angle is rather tight, we will not consider it here. 
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Figure 1: Branching ratios for various decays as a function of £. The solid curve gives the SM 
value. In the presence of a d 222 operator with a sfermion mass of 200 GeV, the long-dashed, 
short-dashed and dot-dashed curves correspond to the cases where each of the two X's equal 0.09, 
0.07 and 0.05 respectively. The thick lines correspond to the experimental bounds. 



a 112- 

Let us first focus on d 222 . In all our subsequent discussions, we take, without any loss of 
generality, both the A's in the product d R to be equal, and the intermediate i-th sneutrino 
mass to be 200 GeV. As is evident from eq. [Ja|, the new physics contribution is proportional to 
A /2 /m~ . Since the products A / 122 A' 1 23 and A 222 A 223 have a stronger experimental upper bound 
than the numbers we need, the only possible solution is for % = 3, i.e., A^Ag^. Similar 
conclusions follow for <if 12 . 

In Figs. (|l]a) and (|1^), we show the effect of a non-zero d 222 on two particular BRs, namely 
those for — > rj'K and 5 ± — > r/K^. Clearly, a resolution of the anomaly is now possible, 
albeit for a A '-dependent range for £. Since the fi p contribution to the decay amplitude tends to 
become too large with increasing A', progressively larger values of £ are required. As for the other 
modes, it is easy to see that the our solutions respect the experimental numbers/constraints. For 
example, with A' = 0.09(0.07) and £ = 0.2(0.3), we expect BR(B° -> t]'K°) = 5(5.5) x 10" 5 , well 
in consonance with observations (Table |I|). Similarly, the BRs for the modes B + — > rjK* + , rj'K* + 
and B° -> r)K°, r]K*°, r)'K*° for A' = 0.09 are predicted to be 1.2(0.6), 0.5(0.3), 0.8(0.4), 0.9(0.4), 
and 0.3(0.2) (xl0~ 5 ) respectively for £ = 0(0.5). In fact, if our explanation be the correct one, 
we would expect to see the decay B — > r]K quite soon, whereas some of the other modes 
may be visible in the upcoming B-factories. 

At this stage, a comment is in order. For Figs.(|I]a, b, c), we have used F = 0.36 and 
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Figure 2: 77&e raizo BR{B ± -»• ri'K ± )/BR(B G -»• vr+A'") /or £ = as a function of the CKM 
parameter 7 . TTie so/id curves represent the SM prediction while the dashed curves are for a d 222 
with each A' = 0.09. In each case, the upper and lower curves are for m s (lGeV) = 150(165) 
Me V respectively. 



m s (lGeV) = 165 MeV (m s (m&) = 118 MeV), values preferred by the SM fit. However, in 
the presence of additional interactions, one may use a different set. As Fig. ^] shows, the depen- 
dence on m s is marginal. On the other hand, a larger value for F would enhance the BRs. For 
example, for F = 0.4, £ = 0.55 and each A' = 0.09, the theoretical BRs for the modes r]'K + , 
rj'K , r]K + , rjK , ir~K + , ir + K°, 7r + 7r~ and 7r + 7r° (last four modes do not have any contribution 
from af 22 ) are 4.9, 7.6, 0.6, 0.7, 2.1, 2.8, 1.1 and 0.9 respectively (all in units of 10~ 5 ). This is 
the maximum value of F that can be used in conjunction with £ = 0.55 since (i) the prediction 
for the 7T~A + mode actually saturates the experimental number, and (ii) the data on B — > nn 
implies that F\V u b/V c b\ < 0.024 (note that semileptonic decays give \V u b/V c b\ = 0.08 ± 0.02). 
Of course, the above does not preclude smaller values for F: with F = 0.33, £ = and each 
A' = 0.09, the theoretical predictions for the abovementioned eight modes are 7.8, 7.6, 0.9, 
0.7, 1.8, 2.9, 1.1 and 0.8 respectively (again, all in units of 10~ 5 ). Anyway, for these numbers, 
particularly for those in the first set, one can easily see that more and more channels get close 
to the discovery limit. 

What about the B — > PV modes? As Fig. [IJ(c) shows, the SM fit requires £ < 0.23. This 
is in conflict with other PV modes such as B ± — * uK^ and B ± —> uix^. The former requires 
either £ < 0.05 or 0.65 < £ < 0.85 while the latter requires 0.45 < £ < 0.85 ||. Interestingly, 
the o?222 operator affects B^ — > (pK^ while the other two decay modes are blind to it. Since 
this additional contribution interferes destructively with the SM amplitude, BR^B^ — > 0A ± ) 
is suppressed leading to a wider allowed range for £ (see Fig |TJc) . For example, with A' = 0.09, 



S 



£ can be as large as 0.8, thus allowing for a common fit to all the three (PV) modes under 
discussion^. d 222 also affects a VV decay modes such as (B — > (f)K*). As this calculation 
involves a few more model dependent parameters, we do not analyse it here. 

Finally, we investigate the consequences for a non-zero rff 12 as opposed to d 22 2- F° r brevity's 
sake, we present graphs (see Fig. ^) only for BR(B ± — > rj'K ). It is interesting to note 
that A' > 0.05 is not admissible for any £ < 1, as the model predictions become significantly 
larger than the observed width. As for B° — > rj'K , the BR is 6.2(4.8) x 10~ 5 for £ = 0.3 and 
A' = 0.025(0.02) (see Table [l]). Indeed, the entire parameter space allowed by B — > rfK is also 
allowed by B° — > rj'K . For such values of A's, BR(B ± — > r/K^ ^ 3x 10~ 6 , and thus well below 




o 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 

0.2 0.4 , 0.6 0.8 1.0 



Figure 3: Branching ratio for B^ — > rj'K ± as a function of £. The solid curve gives the SM 
value. In the presence of a df l2 operator with a sfermion mass of 200 GeV, the long-dashed, 
short-dashed and dot-dashed curves correspond to the cases where each of the two A's equal 
0.025, 0.02 and 0.01 respectively. The thick lines correspond to the experimental bounds. 

the experimental upper limit. Similarly, for the other relevant PP modes B° — > t]K , tt K° 

and B^ —>■ K^ir , the maximum BRs are 0.15, 1.9 and 1.4 (xlO -5 ) respectively. Since, for 

all these decays, the Ifip contribution interferes destructively with the SM one, the resultant 

predictions are considerably suppressed. The best constraints emanate from BR^B^ 1 —>■ K ^) 

which supports 0.03 < £ < 0.8 for the A's used in Fig. [| 

The case for the PV modes is similar. For the decays B + — > i]K* ± ,ri'K* ± ,7r K* ± and 

B° — > r]K*°, r]'K*°, n K*° the $, p operator adds constructively whereas for B^ — > K°p ± , ^K* 

the interference is destructive in nature. The maximum possible BRs for the first six modes, 

for A' = 0.03 and£ = 0(0.5), are 1.2(1.0), 0.75(0.45), 0.4(0.3), 1.7(1.1), 1.0(0.5), 0.8(0.4) (xl0~ 5 ) 

6 Note that the favoured value of £ for the PP and PV modes still continue to be different. While this is not 
a discrepancy, a common £ for both these sets can be accommodated for values of A' slightly larger than that 
we have considered. 
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respectively, smaller than the corresponding experimental numbers. For the last two modes, of 
course, no question of contradiction with experiment arises. 

In short, the modes B ± — > riK*^ and B° — > t]K* are close to the discovery limit whereas 



other modes may have to wait for the next generation B-machines. In a subsequent paper |24 
we will discuss the CP violating effect of these Ifl p operators on all these, and other, modes in 
detail. 



5. Conclusion 

To conclude, we have written down all possible Jji v SUSY contributions to the effective 
Hamiltonian for the B ± — > rj'K decay. We have found that only two new terms, each involving 
two A'-type couplings, can raise the BR to satisfy the experimental number. We have shown 
that though these two terms appear in other nonleptonic decay modes of the B meson, their 
BRs always satisfy the experimental constraints in the whole of the allowed parameter space 
of A', vnv ih and £. Modes like r]K* + , rjK*° are close to their discovery limits. Further, one of 
the new contributions allows larger parameter space in £ for the decay B^ — > (pK^, where the 
other observed modes e.g., B^ — > ouK ± and B^ 1 — > ujtt^ 1 can be fit; this is not possible in the 
SM framework. This leads us to believe that i?-decays and upcoming B-factories may be the 
most promising place to look for new physics beyond the SM. 

We thank Amitava Datta and N.G. Deshpande for illuminating discussions. 
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